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The Dirac equation is solved for triangular and hexagonal graphene quantum dots for different 
boundary conditions in the presence of a perpendicular magnetic field. We analyze the influence 
of the dot size and its geometry on their energy spectrum. A comparison between the results 
obtained for graphene dots with zigzag and armchair edges, as well as for infinite-mass boundary 
condition, is presented and our results show that the type of graphene dot edge and the choice of the 
appropriate boundary conditions have a very important infiuence on the energy spectrum. The single 
particle energy levels are calculated as function of an external perpendicular magnetic field which 
lifts degeneracies. Comparing the energy spectra obtained from the tight-binding approximation 
to those obtained from the continuum Dirac equation approach, we verify that the behavior of the 
energies as function of the dot size or the applied magnetic field are qualitatively similar, but in 
some cases quantitative differences can exist. 
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I. INTRODUCTION 

Since its recent discoverjffl, graphene (a single layer of 
carbon atoms) has been attracting a lot of interest, due 
to its unique band structure, which is gapless and ex- 
hibits an approximately linear dispersion relation at two 
inequivalent points of the reciprocal space (labeled as K 
and K') in the vicinity of the Fermi energy. The linearity 
of the band structure allows one to describe the carriers 
close to the K and K' points in a continuum model, using 
the Dirac equation with massless particles*^. Because of 
the well known Klein tunneling effect in graphene, which 
prevents electrical confinement of electrons, the lateral 
confinement of Dirac carriers is a big challenge in manu- 
facturing graphene-based electronic devices^iEl Different 
suggestions have been made to realize lateral confinement 
of electrons in graphene, e.g. by means of gap engineer- 
ing, provided by a space dependent mass term,!^ or, al- 
ternatively, by combining an external magnetic fieldS or 
a finite mass term^ with an electrostatic potential. On 
the other hand, recent improvements of different fabrica- 
tion techniques made possible cutting and manufacturing 
of single layer graphene flakes, with different shapes and 
size&i2-"— , where such a lateral conflnement naturally oc- 
curs. Using the tight-binding model (TBM), remarkable 
effects have been reported as a consequence of the type 
of the edges and the geometry of these flake^^^HUl {j 
zero-energy states are predicted for triangular graphene 
flakes with zigzag boundaries, ii) for very small flakes a 
gap opens (the energy gap of different graphene flakes 
was recently investigated experimentally and the den- 
sity of states (DOS) strongly depends on the type of the 
edges for any dot geometry, and Hi) the energy levels of 
graphene quantum dots in the presence of a magnetic 
fleld approach the Landau levels with increasing mag- 
netic fleld. 

Recently, analytical results were reported for infinite 




FIG. 1: (Color online) The lattice structure of triangular 
(upper panels) and hexagonal (lower panels) graphene quan- 
tum dots with (a,c) armchair edges and (b,d) zigzag edges. 
a = 0.142 nm is the C-C distance and the primitive lattice 
vectors are denoted by a and b. The atoms of the two sublat- 
tices are represented by blue circles and red dots. The yellow 
region indicates the area of one carbon hexagon. Ng is the 
number of C-atoms in each side of the dot. 



mass boundary conditions for circular disks,'2il for trian- 
gular flakes with armchair edge^^H^nd zigzag edgei^^and 
for square graphene quantum dotJ^. However, it is not 
always clear how the complicated boundary conditions 
describing the zigzag and armchair edges can be invoked 
in the continuum model. Furthermore, the geometry of 
the triangular and hexagonal graphene flakes, make such 
systems harder to be studied by analytical means. One 
has to rely on, either a tight-binding model or a numer- 
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ical solution of coupled differential equations in case of 
the continuum model. 

The continuum model describes very well the low en- 
ergy states in an infinite graphene sheet, but it is not 
clear if this is still the case for small graphene flakes. 
Therefore, it is important to learn if there is a minimum 
size beyond which the continuum model no longer gives 
reliable predictions. Furthermore, because of the large in- 
fluence of the type of edges on the energy spectrum, and 
since it is not always clear which boundary conditions 
should be invoked in the Dirac equation for each possible 
geometry of the flake, a comparison between the results 
obtained with the different possible boundary conditions 
and a link with the TBM is an interesting issue which 
requires a detailed study. 

In this paper, by solving the Dirac equation numeri- 
cally, we present a theoretical study of the energy spec- 
tra of triangular and hexagonal graphene quantum dots, 
where three types of boundary conditions are invoked, 
namely, zigzag, armchair and infinite mass boundary con- 
dition. The infiuence of an external magnetic field, per- 
pendicular to the graphene layer, on the energy spectrum 
of the quantum dots is also analyzed. A comparison be- 
tween the results obtained with the continuum model and 
those obtained from the tight-binding approach will be 
made. 

This paper is organized as follows. In Sec. II we 
present a brief outline of the tight-binding model (TBM) . 
The model based on the Dirac- Weyl equation is presented 
in Sec. Ill and the different boundary conditions are sep- 
arately analyzed in this section. Our numerical results 
are reported in Sec. IV. The summary and conclusions 
of this work are presented in Sec. V. 



II. TIGHT-BINDING MODEL 

The tight-binding Hamiltonian within the nearest 
neighbor approximation is 



4, 



+ h.C.), (1) 



where is the energy of the n-th site, tn^m is the 
hopping energy and c\ (c„) is the creation (annihila- 
tion) operator of the tt electron at site n. Note that, 
for each site n, the summation is taken over all near- 
est neighboring sites m. In the presence of a magnetic 
field, the transfer energy becomes t — > te'^'^*" '", where 
'^n,m — (l/'&o) X''" -A- ■ dl is the Peierls phase, with 
$0 = h/e the magnetic quantum flux and A the vector 
potential. 

Triangular and hexagonal quantum dots with zigzag 
and armchair edges are illustrated in Fig. 1, where the 
vectors a = 0(3/2,^3/2) and b = a(3/2, -^/3/2), with 
a = 0.142 nm the lattice parameter (or the C-C dis- 
tance), are introduced as primitive lattice vectors. In the 
present work, we will consider only the interaction be- 
tween each atom n and its three flrst nearest neighbors. 



In the case of graphene, this interaction has the hopping 
energy t = 2.7 eV. The vector potential corresponding to 
the external magnetic field B = Bz perpendicular to the 
layer is chosen as the Landau gauge A — (0, Bx, 0). With 
this choice of gauge, the Peierls phase for a transition 
between two sites n and m is '^n,m = in the x direc- 
tion and ^n,m = i(2;/3a)<I>c/$o along the ±y direction, 
where <I>c = 3-\/3a^S/2 is the magnetic fiux threading one 
carbon hexagon (the area of one carbon hexagon is shown 
in Fig. [ija) by the yellow region). An external potential 
is represented by a variation in the on-site energies E^, 
and a vacancy or defect can be represented by setting the 
energy of the vacant site to a larger value and the hop- 
ping terms to these atoms as zercpil The Hamiltonian 
H in Eq. Q can be represented in matrix form and the 
eigenvalues and eigenfunctions of a graphene fiake can be 
obtained by diagonalization of the matrix. 

Notice that the hexagonal lattice presented in Fig. 1 is 
not a Bravais lattice, but a combination of two triangular 
lattices composed by atoms labeled as type A (blue) and 
type B (red) . Accordingly, the tight-binding Hamiltonian 
of Eq. ([T]) can be rewritten as 

n n <n,m> 

(2) 

where the operators a\ (an) and (6„) create (annihi- 
late) an electron in site n of lattice A and B, respectively. 



III. CONTINUUM MODEL: 
EQUATION 



DIRAC- WEYL 



Considering an infinite (periodic) graphene sheet and 
after, performing a Fourier transform on the operators 
in Eq. ([T]) and diagonalizing the resulting Hamiltonian 
leads to an energy dispersiorP 



E{\l) = ±t 
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(3) 

The first Brillouin zone in reciprocal space is a hexagon 
with six Dirac points, where only two of them are in- 
equivalent. From the primitive vectors, we can find 
the position of these as X = (27r/3a, 27r/3-\/3a) and 
is:' = (27r/3a, -27r/3y3a). The states near these points 
have approximately a linear dispersion and can be de- 
scribed as massless Dirac fermions by the Hamiltonian 



H = 



Hk 
Hk' 



(4) 



where Hk {Hk') is the Hamiltonian in the K {K') point, 
which are given by 



Hk = VF cr. p, 



H 



K' 



VF a- . p, 



(5a) 
(5b) 
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where a = {ax,(7y) are Pauli matrices and <t* = 
{(Tx,—<Ty) denotes the complex conjugate of the matrix 
(7. In the presence of a magnetic field B perpendicular 
to the graphene layer and using the Landau gauge, one 
can simply rewrite Eq. Q in the following form: 



H = 
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(6) 



where. 



n+ — —ihvT. 



d d 2nB 



(7) 



dx dy $0 
The wave function in real space for the sublattice A is 

and for sublattice B it is given by. 



(8b) 



The Hamiltonian of Eq. ([6| acts on the four-component 
wave function = [ipA,ipB,fA',^B']'^ , which leads to 
the four coupled first-order differential equations: 
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In the above equations, we used the following dimension- 
less units: x' = xj^fS, y' = y/v^, /3 = 27rSS'/$o = 
27r$/$o- e = EIEq, with £'0 = Hvp/VS, where S (x L'^ 
is the area of the dot with L being the length of the 
side of the dot. In this paper, we solve Eqs. ^ nu- 
merically, using the finite elements method, for the tri- 
angular and hexagonal graphene fiakes shown in Fig. 1, 
considering zigzag, armchair and infinite mass bound- 
ary conditions. The numerical calculations are performed 
by using the standard finite element package COMSOL 
Multiphysicd^H^ which discretizes the two-dimensional 
fiake in a finite-sized mesh and allows the implementation 
of the appropriate boundary conditions. The way the 
boundary conditions are implemented in the continuum 
model is the subject of the following three subsections. 



A. Zigzag boundary conditions 

The geometry of the hexagonal and triangular 
graphene quantum dots with zigzag edges are illustrated 



VB' = (iV3/2 + l/2)v3A' 



fB = (j^/2 - 1/2)^A 
VB' = (-i\/3/2 + l/2)ipA' 




FIG. 2: The infinite-mass boundary conditions implemented 
on the edges of a triangular dot. h-i,h2,n-i are the outward 
unit vectors at each edge of the dot. 
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FIG. 3: Energy levels of hexagonal (a, c, e) and triangular 
(b, d, f) graphene quantum dots with zigzag (a, b), armchair 
(c, d) edges and infinite mass boundary condition (e, f) as 
function of the square root of the dot area S in the absence 
of a magnetic field. 



in Figs. l(b,d). The length of one side of the hexago- 
nal and triangular dots, respectively, are given by L = 
V^iNs - l/3)a and L = V^{Ns + l)a, with being the 
number of atoms in each side of the dot and a = 0.142 nm 
is the C-C distance. The total number of C-atoms in the 
triangular dot is = [(A^s + 2)^-3] and N = 6N^ for the 
hexagonal dot. The zigzag- type boundary condition was 
previously studied by Akhmerov et al?^ , who presented a 
model which is generically applicable to any honeycomb 
lattice. For a graphene dot with zigzag edges and if the 
last atoms at the boundary are from sublattice A (blue 
circles in Fig. 1), the boundary conditions are given by 
LpA = ^A' = 0, whereas Lpg and LpB' are not determined, 
and similarly, when the zigzag edges are terminated by 
the B atoms (red dots in Fig. 1), ips — fB' = 0, while 
ipA and ifA' are not determined. 
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FIG. 4: Energy levels of a zigzag hexagonal graphene dot as 
function of the eigenvalue index obtained by (a) the TBM and 
(b) the continuum model, for three different sizes of the dot 
with Ns = 10, 20, 30, having respectively surface area S = 
14.68, 60.78, 138.32 nm^. The inset in panel (a) shows the 
energy gap Eg as function of A'^a obtained, by the TBM. 



B. Armchair boundary conditions 
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FIG. 5: Energy levels of an armchair hexagonal graphene dot 
as function of the eigenvalue index obtained by (a) the TBM 
and (b) the continuum model for three different sizes of the 
dot with Ns = 20, 40, 60 having respectively surface area 
5=41.07, 176.23, 405.68 nm^. The inset in panel (a) shows 
the energy gap obtained from both TBM (black squares) and 
continuum model (green circles). The inset in panel (b) shows 
the lowest electron energy levels as function of Ns for both 
TBM (blue solid curves) and continuum model (red dashed 
curves). 



The geometry of a hexagonal and triangular graphene 
quantum dot with armchair edges is illustrated in Figs. 
l(a,c) . Here, the length of one of the edges of the hexagon 
dot is L — {3Ns — 4) a/2 and for the triangular dot is 
L = 3Nsa/2. For an armchair hexagonal graphene dot 
the total number of C-atoms is iV = [9Ns{Nj2 - 1) + 6] 
and for the triangular dot is given by = {Ns+2)3Ns/4:. 
Note that in the case of armchair boundaries the number 
of C-atoms in each side is an even number (see Figs. 
[l];a,c)). 

From Figs. l(a,c), we notice that the edge atoms con- 
sist of a line of A-B dimers, where the wave function 
should be zero. From Eqs. 



ary conditions becom^^ZI 



and (8b), these bound- 



- _A^'-K)-r 



(10a) 
(10b) 



where r is taken at the position of the edge. Notice 
that these armchair boundary conditions mix the wave 
functions of the K and K' points. 



C. Infinite-mass boundary condition 

A mass-related potential energy V{x, y) can be coupled 
to the Hamiltonian via the cr, Pauli matrix, 



H = vp a. p 



'y{x,y), 



(11) 



where the parameter t — ±1 distinguishes the two K 
and K' valleys. It is straightforwardly verified that the 
presence of a mass term in the Hamiltonian of Eq. (11 1 



induces a gap in the energy spectrum of graphene. How- 
ever, if the mass-related potential V{x,y) is defined as 
zero inside the dot and infinity at its edge, the Klein 
tunneling effect at the interface between the internal and 
external regions of the dot can be avoided and, con- 
sequently, the charge carriers will be confined. This 
infinite-mass boundary condition can be introduced in 
the Dirac equation by defining (pB{x,y)/(pA{x,y) = ie^'^ 
and (pBr{x,y)/(pA'ix,y) = —ie^'^ (which, respectively, 
correspond to the if-point and the ii''-point wavespinors) 
at the boundary, where </> is the angle between the out- 
ward unit vector at the edges and the x-axis.l^Due to its 
simplicity, this type of boundary condition has be en us ed 
in the study of circular graphene dot^Sil and ringd^HES! in 
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FIG. 6: Electron probability densities corresponding to the 
two lowest energy levels of hexagonal graphene flakes, ob- 
tained by the continuum model, for (a,b) zigzag {Ns = 20) 
and (c,d) armchair {Ns = 40) edges. 
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FIG. 7: Energy levels of a zigzag triangular graphene dot as 
function of the eigenvalue index obtained by (a) the TBM and 
(b) the continuum model for three different sizes of the dot 
with Ns = 12, 24, 40 having respectively surface area 5=4.42, 
16.37, 44.03 nm^. The inset in panel (a) shows the energy 
gap obtained from both TBM (black squares) and contiimum 
model (green circles). 
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FIG. 8: Energy levels of an armchair triangular graphene dot 
as function of the eigenvalue index obtained by (a) the TBM 
and (b) the continuum model for three different sizes of the 
dot with Ns = 20, 40, 60 having respectively surface area 
S=7.85, 31.43, 70.72 nm^. The inset in panel (a) shows the 
energy gap obtained from both TBM (black squares) and con- 
tinuum model (green circles). The inset in panel (b) shows 
the lowest electron energy levels as function of Ns for both 
TBM (blue solid curves) and continuum model (red dashed 
curves). 
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FIG. 9: Electron probability densities corresponding to the 
lowest energy levels of the triangular graphene flakes, ob- 
tained by the continuum model, for (a,b) zigzag {Ns = 20) 
and (c,d) armchair {Ns = 40) edges. 
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FIG. 10: Electron densities for the first energy level of the 
triangular and hexagonal graphene flakes (using TBM) with 
Ns — 10 and zigzag edges. Left panels show the total electron 
density j*!*]^ and the right panels present the electron densities 
associated with A and B sublattices. The gray dots are the 
positions of C-atoms. 



the presence of a perpendicularly magnetic field, where 
analytical solutions can be found. For the hexagonal and 
triangular geometries the angle (f> has a fixed value at 
each side of the dot that simplifies the boundary condi- 
tions to ifB = ctipA (for the K valley) and ipB' = —ctipA' 
(for the K' valley) where a = ie**^ is a complex num- 
ber. The infinite-mass boundary conditions are shown 
explicitly in Fig. |2]for a triangular dot. 



IV. NUMERICAL RESULTS 

A. Zero magnetic field 

The energy levels of hexagonal (upper panels) and 
triangular (lower panels) graphene flakes, as calculated 
within the continuum model, are shown in Fig. [3]as func- 
tion of the square root of the dot area. The results are 
shown for zigzag (a,b), armchair (c,d) and infinite-mass 
(e,f ) boundary conditions and are qualitatively and quan- 



FIG. 11: The same as Fig. [Tojbut for the dots with armchair 
boundaries and Ns =20. 



titatively very different. As the dot area increases, the en- 
ergy levels tend to a gapless spectrum, which is expected, 
since the energy spectrum of an infinite graphene sheet 
does not exhibit a gap. A peculiar spectrum is observed 
for zigzag triangular dots (Fig. [sj^b)): zero energy states 
are found for all sizes of such a dot. These zero energy 
states are separated from the remaining positive and neg- 
ative energy states by an energy gap which decreases as 
the dot becomes larger. The presence of such zero energy 
states in triangular and trapezoidal graphene flakes have 
been previously reported in the literaturd^^"^, where the 
TBM was applied. In the case of zigzag triangular dots, it 
has been shown analyticalljP^ that the equation = 
for the TBM Hamiltonian in Eq. Q leads to iV^ - 1 
linearly independent states, namely, Ns — 1 degenerate 
states with E = 0, for any number of C-atoms in 
one of the edges of the flake. Thus, Fig. |3jb) demon- 
strates that the existence of zero energy states, which is 
observed in the TBM, is qualitatively captured by the 
approximations of the continuum model as well. The 
results in Fig. |3] also show that the energy levels for a 
dot with armchair and infinite-mass boundary conditions 
are qualitative more similar to each other than the spec- 
tra for zigzag edges, where carriers are predominantly 
confined at the edge of the dot. In fact, for the trian- 
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FIG. 12: Energy levels of (a,c) triangular and (b,d) hexago- 
nal graphene dots with armchair boundary as function of the 
magnetic flux threading one carbon hexagon <l>c- The results 
in panels (a,b) are obtained using the continuum model while 
panels (c,d) display the TBM results. The quantum dots have 
an area S such that \/S — 10 nm. 

gular geometry, the infinite mass boundary condition de- 
scribes very well the armchair states, specially for lower 
energy states. However, for the hexagonal geometry, the 
results for armchair and infinite mass boundary condi- 
tions are only qualitatively similar where, the hexagonal 
dots with infinite mass boundary condition exhibit more 
energy states in comparison with the armchair case. 

Notice that the energy spectra shown in Fig. [3] ex- 
hibits degenerate states. These degeneracies, which will 
be evidenced in the following figures, where we plot the 
energy spectra as a function of the eigenvalue index, are 
related to the symmetries of the triangular and hexag- 
onal dots, as we will explain in further detail later on, 
when we discuss about the electron probability densities. 

A comparison between the energy spectra obtained by 
means of the TBM (a) and the Dirac equation (b) for 
zigzag hexagonal dots is shown in Fig. |4] for three sizes 
of the dot, defined by the number of C-atoms in each 
side of the hexagon Ng . The energies Ei are plotted as a 
function of the eigenvalue index i. Although the results 
are quantitatively different, they are qualitatively simi- 
lar, e.g. as the size of the dot increases, they start to 
exhibit an almost fiat energy spectrum as a function of 
the eigenvalue index around the Dirac point. Such a fiat 
spectrum leads to a peak in the DOS close to the Dirac 
point, which was recently reported in the literature^ for 
graphene dots with zigzag edges within the TBM. The 
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FIG. 13: The same as Fig. |12| but for dots with zigzag bound- 
aries and the dots have an area S such that \^ = 5 nm. 



curves for Ng — 30 obtained by the TBM and continuum 
models are very similar, except for the fact that many 
more states are found in the latter, whereas the discrete 
character of the spectrum in the former is much more 
clear. For smaller dots, the agreement between these two 
models becomes clearly worse. For instance, an energy 
gap Eg is found for very small hexagons (i.e. Ng < 10) 
within TBM, whereas in the case of the continuum model 
such a gap is extremely small. As a consequence, the con- 
tinuum model overestimates the DOS at £^ = as the dot 
size decreases, since it exhibits a plateau in the energy 
as a function of the eigenstate index in the vicinity of 
E = even for smaller Ng, where TBM results show 
a gap in the energy spectrum. Notice that the E = 
states in zigzag dots are edge states, so that the num- 
ber of zero-energy states depends on the number of edge 
atoms in the TBM and, similarly, to the number of mesh 
elements at the edge in the continuum model. Therefore, 
in the continuum model for E — 0, the finite elements 
problem is ill-defined, where the constructed matrix of 
the finite mesh elements in this case is singular (zero in- 
verse), leading to spurious solutions around E = 0. As 
the size of the dot increases, the gap in the TBM results 
quickly reduces to zero and a zero energy level for the 
hexagonal fiakes with zigzag edges appears.— In the in- 
set of Fig. Qa) , the energy gap values obtained by TBM 
are shown as function of Ng. These results can be fit- 
ted to Eg = a{l/Ngy (blue solid curve in the inset of 
Fig. ^a.)), where a = 94.6 eV and 7 = 3.23 are fitting 
parameters. 
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FIG. 14: The same as Fig. |12[ but for dots with infinite-mass 
boundary conditions. The dots have an area S such that 
\/S = 5 nm. In panels (c,d), the infinite-mass boundary con- 
dition is applied within the TBM model, where we imposed a 
+10(-10) eV on-site potential for sublattice A(B) around the 
dot geometry. 



The energy states of armchair hexagonal dots are 
shown as a function of the eigenvalue index in Fig. [5] 
within the TBM approach (Fig. 5ja)) and the Dirac- 
Weyl equations (Fig. [5](b)), for three different sizes of 
the dot. The energy spectrum in both cases approach the 
prolonged S-shape curve predicted by Ezaw#31 as the size 
of the dot increases and the spectrum exhibits an energy 
gap Eg at the Dirac point. The energy gap as function 
of Ng is shown in the inset of Fig. [5](a) which decreases 
rapidly as the size of the dot increases. Our numerical 
results can be fitted to Eg = a/Ng with a = 8.5 eV for 
the TBM (blue solid curve) and a = 13 eV for the con- 
tinuum model (red dashed curve) results. Notice that 
Eg obtained from the continuum model is larger than 
the one from the TBM results in particular for small Ng 
and both curves can not be made to coincide by a simple 
shift in Ng. This is clearly a consequence of the increased 
importance of corrections to the linear spectrum used in 
the continuum model for small sizes of the system. The 
inset of Fig. [sjb) shows the five lowest electron states for 
both TBM (blue solid curves) and the continuum model 
(dashed red curves). Our results show that the contin- 
uum model overestimates the energy values also for the 
upper energy levels in comparison with the TBM energy 



levels. In fact, the energy dispersion in the continuum 
model is given by a linear curve, which coincides with 
the TBM energy spectrum for low energies, but as the 
energy goes further away from E = 0, this linear dis- 
persion overestimates the energy as compared to the real 
band structure of graphene, which starts to bend down 
from the linear spectrum as the energy increases. This 
emphasizes once again the importance of the higher or- 
der corrections to the linear dispersion, especially for high 
energy states and smaller dot sizes. 

Figure [6] shows the probability density (using the con- 
tinuum model) corresponding to the first two energy lev- 
els of hexagonal flakes. The probability density for the 
zigzag case with Ng — 20 is presented in panels (a) and 
(b), respectively, for E = and E = 0.01 eV. The re- 
sults clearly demonstrate that the zero energy states in 
the zigzag case are due to edge effects and, accordingly, 
are confined at the edges while the carriers confine to- 
wards the center of the flake with increasing energy (see 
Fig. |6]^b)). The probability densities of the armchair 
edged graphene flake with Ng = 40 are very different as 
seen in Figs. |6jc,d) for the lowest degenerate states with 
E — 0.16 eV. The electron wavefunction is spread out 
over the whole sample, but different from the usual quan- 
tum dots with parabolic energy-momentum spectrum, it 
has a local minimum in the center of the dot. Note that 
Fig. |6jc) has only two-fold symmetry while Fig. [6](d) is 
six-fold symmetric. Both densities are zero in the center, 
while Fig. [gJ^c) has two extra zeros at the sides along 
y — 0. These results are comparable to the TBM results 
obtained in ReflTZl 

The energy spectrum for triangular dots with zigzag 
edges, obtained by the TBM and the Dirac- Weyl equa- 
tion are shown as a function of the eigenvalue index in 
Figs. |7][a) andjrf^b), respectively. Notice that both en- 
ergy spectra exhibit zero energy states. As we mentioned 
before, the number of degenerate states with zero energy 
is a well defined quantity in the tight-binding approach, 
namely, Ng — 1, where Ng is the number of C-atoms in 
one side of the triangl^^. On the other hand, the re- 
sult in Fig. [rj^b) for the continuum model exhibits many 
more zero energy states. Therefore, while the contin- 
uum model captures qualitatively the existence of zero 
energy states, it does not provide the appropriate num- 
ber of degenerate states as calculated by the TBM. These 
zero energy levels are related to the edge states of zigzag 
graphene flakes^** The energy gap (between the zero 
energy level and the flrst non-zero eigenvalue) is shown 
in the inset of Fig. |7|^a) as function of the size of the 
dot, where Eg obtained by both models are comparable 
and the difference between the TBM (red dashed curve) 
and continuum (blue solid curve) results tends to zero 
for large graphene flakes. These results can be fitted to 
Eg = a/Ng with a = 15.75 eV for the TBM gap and 
a — 18.9 eV for the continuum model. 

The energy spectra of triangular dots with armchair 
edges obtained by the TBM and the continuum model 
are shown in Fig. [8] No zero energy states are found 
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FIG. 15: The same as Fig. 3 but in the presence of an external 
magnetic field of B = 50 T. The inset shows the energy gap 
as function of the magnetic flux obtained by the TBM for two 
values of A^s. The triangle and circle symbols display Eq. (13) 
which is fitted to the numerical results. 
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FIG. 16: The same as Fig. 4 but in the presence of an external 
magnetic field of _B = 50 T. The inset shows the energy gap 
as function of the magnetic flux obtained by the TBM (solid 
curves) and continuum model (dashed curves) for two values 
of A^s- The triangle and circle symbols display Eq. (13) which 
is fitted to the numerical results. 



and the energy gap at the Dirac point for both models 
is comparable. The gap can be fitted to Eg = a/Ng 
{a ^ 21.9 eV for TBM and a = 25.9 eV for the contin- 
uum model) as shown respectively by the blue solid and 
dashed red curves in the inset of Fig. [sj^a). The lowest 
electron energy levels, obtained by the TBM (blue solid 
curves) and the continuum model (red dashed curves), 
are shown in the inset of Fig. [sf^b) as function of Ns- 
The results show a larger difference between the TBM 
and continuum energy values for the upper energy levels 
(e.g. \Ef-E?\<\E^-E§\). 

Notice that the energy gaps found for all the systems 
that we investigated were fitted to Eg = a/N^ for dif- 
ferent values of a, except for the case of zigzag hexag- 
onal dots, where the gap is fitted to Eg = a/N] . with 
7 = 3.23. This is a consequence of the fact that the 
corners of the zigzag hexagonal dot structure are not ter- 
minated by a single atom, as in the case of zigzag tri- 
angular dots, but by a pair of C-atoms corresponding 
to two different sublattices, forming a A-B dimer (see 
Fig. 1). These A-B dimers are responsible for a vanish- 
ing wave function in the corners of the zigzag hexagonal 
dots, as observed in Fig. |6] As explained in Sec. Ill A, 
the zigzag boundary condition for each side of the dot is 
implemented in the Dirac- Weyl equations by setting to 
zero the component of the pseudo-spinor corresponding 
to the sublattice that forms that side. As the sublattice 
types of adjacent sides of a zigzag hexagonal dot are dif- 
ferent, connected by the A-B dimers in the corners, the 
whole wave function must vanish at these corners, since 
these points are composed of both A and B sublattices. 
The vanishing wave function at the corners reduce the ef- 
fective confinement area and, consequently, increases the 
energy gap, especially for smaller dots, where the infiu- 
ence of the corners is more significant. As the size of the 
dot increases, the role of the corners in the energy gap 
becomes less important and is eventually suppressed by 
the influence of the zigzag edges, leading to the zero en- 
ergy states that form the plateau in Fig. |4] explaining 
the faster decay of the energy gap (7 = 3.23) in zigzag 
hexagonal dots, as compared to the other cases (7 — 1). 

The probability density corresponding to the first two 
energy levels of triangular graphene flakes, obtained by 
the continuum model, is shown in Fig. [9] The proba- 
bility density for the zigzag edged dot with — 20 is 
presented in panels (a) and (b), respectively, for = 
and for the first non-zero eigenvalue (i.e. E = 0.92 eV). 
For the degenerate zero energy states the carriers are con- 
fined at the edges of the triangular fiake which is typical 
for zigzag boundaries. States corresponding to large en- 
ergy values are confined in the center of the triangle (Fig. 
|9jb)). For armchair triangular fiakes, as in the hexago- 
nal case, the electron state is spread out over the whole 
fiake (Figs. |9jc,d) display the different probability den- 
sities for Ns = 40 corresponding to the first degenerate 
eigenvalues with E = 0.32 eV). Both wavefunctions have 
three-fold symmetry and the inner part is even six-fold 
symmetric. Note that the electron density in Fig. l9|d) is 
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Figs. 
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FIG. 17: The same as Fig. 5 but in the presence of an external 
magnetic field of _B = 50 T. The inset shows the energy gap 
as function of the magnetic flux obtained by the TBM (solid 
curves) and continuum model (dashed curves) for two values 
of Ns- The triangle and circle symbols display Eq. (13) which 
is fitted to the numerical results. 



zero at the three corners and in the center of the triangle 
which is different from Fig. ^c) where zero's are found 
at the corners of the inner hexagon and at the center of 
the sides of this hexagon. 

The TBM electron densities of the zigzag graphene 
dots with Ns = 10 is shown in Fig. [TO] for the first 
energy level of the triangular and hexagonal graphene 
flake. Left panels present the total electron density \'^'^\ 
and the electron densities associated with A and B sub- 
lattices (|pA,Bp) are shown in the right panels. We found 
that the wavefunctions of the two-fold degenerate states 
are related to each other by a 60° rotation. The sum 
of the densities of the degenerate states results in a six- 
fold (three-fold) symmetric wavefunction for the hexag- 
onal (triangular) flakes. As seen in Fig. flOl the total 



As seen in Fig. 
electron density is related to the densities of A and B 
sublattices by l^fp — \pa\^ + \pb\'^ ■ Figure 11 describes 
the density distributions of the lowest energy levels for 
armchair graphene flakes. For the armchair hexagonal 
dots the electron densities corresponding to the A and 
B sublattices (right panels) can be transformed to each 
other by a 180° rotation whereas the density of the tri- 
angular wavespinors can not be linked to each other by 
a rotational transformation. 



B. Magnetic field dependence 

The dependence of the energy levels of triangular (a) 
and hexagonal (b) graphene flakes on the magnetic flux 



12 and 13 respectively for flakes with armchair and 
zigzag edges. The results in panels (a),(b) are obtained 
using the continuum model and the results in panels 
(c),(d) show the TBM energy spectrum. Sc = (3%/3/2)a2 
is the area of a carbon hexagon which is indicated by the 
yellow region in Fig. ^a). The results are obtained for 
dots with an area of 5 = 100 nm^ and S = 25 nm^ re- 
spectively for armchair and zigzag edges. The continuum 
and TBM results are qualitatively similar to each other 
in the sense that as the magnetic flux increases, the en- 
ergy levels converge to the Landau levels of a graphene 
sheet En (see red solid curves), which are given by 



En = sgn{n) 



(12) 



where. Is = \/ h/eB is the magnetic length and n is an 
integer. The interplay between the quantum dot con- 
flnement and the magnetic fleld conflnement is respon- 
sible for the appearance of a series of (anti)-crossings 
in the energy spectrum. As explained earlier, armchair 
graphene dots do not exhibit zero energy states for B — 0. 
However, as the magnetic field increases, some of the ex- 
cited energy levels approach the zero energy Landau level 
n = in both armchair and zigzag graphene fiakes, which 
naturally produces (anti)-crossings between the excited 
states. Lifting the degeneracy of the energy levels by the 
magnetic field results in a closing of the energy gap with 
increasing magnetic field. Notice that the zero energy 
states of zigzag triangular dots (Fig. [Tsf^a)) are not af- 
fected by the magnetic field because they are strongly 
confined at the edges of the dot. All these features are 
qualitatively similar to th ose ob tained by the TBM (see 
the lower panels in Figs. 12|13[ ). In the case of hexago- 
nal zigzag graphene dots (see Fig. |13[ b)), the continuum 
model exhibits a plethora of additional lines as compared 
to the well known energy spectrum obtained by the TBM 
(compare Figs. 13 b) and [l3|c)). 

For the infinite-mass boundary condition, the energy 
spectrum of triangular (a) and hexagonal (b) dots as a 
function of the magnetic field is shown in Fig. [14] for the 
dot with area 5 = 25 nni'^. The energy spectrum in this 
case differs from both obtained for zigzag and armchair 
boundary conditions. The spectra exhibit no zero energy 
state at i? = and show crossings and anti-crossings be- 
tween the higher energy levels which resemble the TBM 
results (see Figs. [l4[c),(d) respectively for triangular and 
hexagonal dots). In the TBM model, the infinite-mass 
boundary conditions can be realized as a graphene dot 
structure surrounded by an infinite mass media, where 
we applied a staggered potential (i.e. +10(-10) eV on-site 
potential for sublattice A(B)) around the dot geometry. 

The energy levels obtained by the TBM (a) and the 
continuum model (b) for hexagonal graphene flakes un- 
der a B = 50 T external magnetic field are shown in 
Figs. [T5| and [16] for zigzag and armchair edges, respec- 
tively, as function of the eigenvalue index. The energy 
spectra of such systems in the absence of magnetic field. 
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which are shown in Figs. |4] and 5] are composed of a 
series of degenerate states for |£' > 0. The magnetic 
field hfts the degeneracy of such states and reduces the 
gap between the states. The energy gap as function of 
the magnetic flux through a single carbon hexagon <i>c is 
shown in the insets of Fig. [isj^a) and Fig. 16 a) respec- 
tively for zigzag (with Ng = 10, 20) and armchair (with 
Ng = 20, 40) hexagonal dots. These results can be fitted 
to 

^<,($e/$o) = ^" + i^g($c/*o) + i;g'($c/$o)' (13) 

where E^^^ (eV) are the fitting parameters. In the inset 
of Figs. ^Sja) and[l6ja); the fitted results are shown by 
symbols. The fitting parameters for the TBM results of 
a zigzag hexagonal dot with Ng — 10 (for the range of 
< $c/*o < 0.17 ) are E'^-'^^^ ^ (0.12,-0.91,1.36) eV 
(see triangles in the inset of Fig. 
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a)) and E°'^' 



(0.86,-26,210) eV, E°''^''^ = (0.88,-12.5,46.5) eV are 
the fitting parameters of an armchair hexagonal dot with 
Ng = 20 respectively for TBM and continuum results 
(triangles in the inset of Fig. [T6[a)). The fittings are 
done for the range of < $c/*o < 0.06 and < ^c/^o < 
0.13 respectively for TBM and continuum results. 

For the zigzag case and for Ng ~ 20, the energy gap is 
already negligible, whereas for Ng — 10, the Eg « 0.12 
eV gap at -B = decays as the magnetic fiux increases 
and approach zero in the limit of large magnetic fiux (i.e. 
$c/$o > 0.2). Due to the lifting of the degeneracies, the 
energy spectrum of an armchair hexagonal dot exhibits 
an almost linear behavior around £^ = as function of 
eigenvalue index where, both TBM and continuum mod- 
els approximately display the same slope for the linear 
regime. 

For triangular graphene fiakes under a _B = 50 T (i.e. 
^c/^o — 0.0063) magnetic field, the energy spectra ob- 
tained by the TBM (a) and the continuum model (b) are 
shown in Fig. [17] considering zigzag edges, and Fig. [TS] 
considering armchair edges. As mentioned earlier, the 
presence of a magnetic field does not affect the E = 
edge states in the triangular zigzag fiakes, but lifts the 
degeneracy of the E states. The energy gap Eg 
around E' = of triangular fiakes is shown as function of 



magnetic fiux $c in the insets of Fig. [T7|a) and Fig. 18 'a) 
respectively for zigzag (with Ng = 12, 24) and armchair 
(with Ng = 20, 40) edges (circle and triangle symbols 
present the fitted results). E°'^'^ = (1.12, -1.32, -0.028) 
and Eg'^''^ — (1.5,-1.77,0.4) are the fitting parameters 
of a zigzag triangular dot with Ng = 12 respectively for 
TBM and continuum results (see inset of Fig. [iT^a)). 
The fitting parameters for an armchair dot with Ng = 20 
(see inset of Fig. [Tst^a) ) obtained by TBM and continuum 
models are respectively E^'^'^ = (1.02,-3.87,3.83) and 
i;0>i.2 ^ (1 12, -2.41, 1.2). In both zigzag (with Ng = 12 
and armchair {Ng = 20) triangular dots the fittings are 
done for the range of < ^c/^o ^ 0.2). In contrast 
with hexagonal dots the energy gap of the triangular dots 
reduces smoothly (i.e. almost linearly) with increasing 
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FIG. 18: The same as Fig. 6 but in the presence of an exter- 
nal magnetic field B = 50 T. The inset in panel (a) shows the 
energy gap Eg, obtained by the TBM (solid curves) and con- 
tinuum model (dashed curves), as function of the magnetic 
flux through one carbon ring $c for TVs = 20 and A''^, = 40. 
The triangle and circle symbols display Eq. (13) which is fit- 
ted to the numerical results. In the inset of panel (b) Eg is 
shown as function Ns for both TBM (black squares) and con- 
tinuum models (green circles) in the presence of an external 
magnetic field B = 50 T. 



the magnetic fiux. Therefore the energy gap is weakly 
affected by a low magnetic field in triangular graphene 
dots. In the inset of Fig. [iSjb) the energy gap is shown 
as function of Ng. As in the case of zero magnetic field 
Eg can be fitted to Eg = a /Ng as function of Ng (see blue 
solid and red dashed curves in Fig. 18 |b)). The fitting 
parameter for S = 50 T are a « 21.87 eV for TBM and 
a « 25.9 eV for the continuum model which is almost the 



same as for zero magnetic field (see Fig. [s]), i.e. because 
of the linear magnetic field dependence of the energy gap 
for low magnetic field it does not affect significantly the 
dependence of the energy gap on the size of the armchair 
triangular graphene dot. 

As a matter of fact, tuning the energy gap by adjusting 
the external magnetic field is more useful for smaller sizes 
of the dot, since the energy gap decays to zero as the size 
of the dot increases. On the other hand, due to the small 
size of the dots considered in Figs. [15]- [18] we need large 
magnetic field values (e.g. B = 50 T) in order to see 
its effect on the energy spectrum. Nevertheless, as the 
influence of the magnetic field scales with the magnetic 
fiux through the dot area, similar results will be obtained 
for lower magnetic fields in case of a larger graphene dot. 
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V. SUMMARY AND CONCLUDING REMARKS 

We have presented a theoretical study of triangular 
and hexagonal graphene quantum dots, using the two 
well known models of graphene: the tight-binding model 
and the continuum model. For the continuum model, 
the Dirac-Weyl equations are solved numerically, consid- 
ering armchair, zigzag and infinite mass boundary condi- 
tions. A comparison between the results obtained from 
the TBM and the Dirac-Weyl equations show the impor- 
tance of boundary conditions in finite size graphene sys- 
tems, which affects their energy spectra. The results ob- 
tained by the TBM for graphene flakes are only qualita- 
tively similar to the results from the Dirac-Weyl equation 
for such systems considering zigzag and armchair bound- 
ary conditions, which shows that energy values obtained 
from the continuum model for small graphene dots may 
not always be quantitatively reliable. 

More specifically, for zigzag hexagonal and triangular 
dots, the DOS at i? = in the absence of a magnetic field 
is overestimated in the continuum approach. Similarly, 
the continuum model also overestimates the energy gap 
around i? = in the armchair case for both geometries. 
A good agreement between both models is only observed 
for very large dots, as expected, and such agreement is 
always better for the triangular case, as compared to the 
hexagonal case. The energy spectrum obtained using the 
continuum model with infinite mass boundary condition 
for hexagonal graphene flakes do not exhibit the same 
properties as the results obtained with the armchair or 
zigzag boundaries (in both TBM and continuum mod- 
els), which shows that this type of boundary condition 
may not give a good description of finite size hexagonal 



graphene flakes. On the other hand, for the triangular 
case, the results from the continuum model with infinite 
mass boundary conditions describe very well the case of 
triangular dots with armchair edges. 

In the presence of an external magnetic field, the 
energy levels obtained by the continuum model with 
zigzag and armchair boundary conditions converge to 
the Landau levels of graphene as the magnetic field 
increases, as observed in the TBM. However, many ad- 
ditional artifact states appear in the continuum model, 
which do not match with any TBM result and do not 
approach any Landau level. Besides, the influence of an 
external magnetic field on the gap in the energy spectra 
of graphene fiakes is particulary different for triangular 
and hexagonal dots. The energy gap of the hexagonal 
fiakes (with < 10) reduces quickly with increasing 
the magnetic fiux, whereas the gap of the triangular 
fiakes decreases smoothly as the magnetic flux increases. 
This feature is observed in both TBM and continuum 
model, and suggests that the energy gaps of hexagonal 
flakes are more easily controllable by an applied ex- 
ternal field, as compared to the triangular graphene dots. 
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